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ABSTRACT 

We present results from deep X-ray stacking of >4000 high redshift galaxies from z ~1 to 8 using 
the 4 Ms Chandra Deep Field South (CDF-S) data, the deepest X-ray survey of the extragalactic sky 
to date. The galaxy samples were selected using the Lyman break technique based primarily on recent 
HST ACS and WFC3 observations. Based on such high specific star formation rates (sSFRs): log 
SFR/M* > —8.7, we expect that the observed properties of these LBGs are dominated by young stellar 
populations. The X-ray emission in LBGs, eliminating individually detected X-ray sources (potential 
AGN), is expected to be powered by X-ray binaries and hot gas. We find, for the first time, evidence 
of evolution in the X-ray/SFR relation. Based on X-ray stacking analyses for z < 4 LBGs (covering 
~ 90% of the Universe's history), we find that the 2-10 keV X-ray luminosity evolves weakly with 
redshift (z) and SFR as: log Lx = 0.93 log(l+z)-|-0.65 logSFR+39.80. By comparing our observations 
with sophisticated X-ray binary population synthesis models, we interpret that the redshift evolution 
of Lx/SFR is driven by metallicity evolution in HMXBs, likely the dominant population in these high 
sSFR galaxies. We also compare these models with our observations of X-ray luminosity density (total 
2-10 keV luminosity per Mpc'^) and find excellent agreement. While there are no significant stacked 
detections at z > 5, we use our upper limits from 5 < z < 8 LBGs to constrain the SMBH accretion 
history of the Universe around the epoch of reionization. 



1. INTRODUCTION 

Recently, Chandra completed the deepest X-ray survey 
to date: 4 Ms in th e Chandra Dee p Field-South (CDF- 
S: iXue et all [2011 hereafter IxTTh . The last leap for- 
ward in X-ray survey depth occurred approximately a 
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decade ago as deep surveys moved fr om 1 Ms to 2 Ms 
coverage {e.g., [Alexander et al.l [ 2003l) . Such deep X-ray 
surveys are crucial to studying the high-energy emission 
from distant star-forming galaxies, which overall are an 
extremely X-ray faint population compared to the more 
luminous active galactic nuclei (AGN) population that 
dominates the nurn ber counts at bright fluxes {e.g., see 
iLehmer et al.l[20T^ . Notably, since the profusion of pa- 
pers several years ago on high redshi ft (z > 2) X-ray 
studi es of star-forming galaxies (i.e.. iRe ddv fc Steidell 
2004: iLaird eFall [20051 [200l ILehmer et al . 2005b) there 
have been significant improvements in the identification 
of new high redshift galaxy samples, selected using the 
Lyman break technique, from z = 1.5—8 jB eckwith et al. 
[2001 iBou wens et all [2007!. I2008L [2OI0I: iReddv fc Steidel 
2009; Oesc h et al.| [20l5b: H athi et al.ll2010l:lBunker et al. 
I2OIQ) . Thus, it is timely to study the evolution of X-ray 
properties in star-forming galaxies with these updated 
data sets. 

The Lyman break technique efficiently selects dis- 
tant star-forming galaxies based on a strong spec- 
tral break at 912 A, caused by Lyman series absorp- 
tion of neutral hydrogen in the intergalactic medium; 
such galax ies are referred to as Lyman break galax- 
ies (L BGs: [steidel fc Hamilton|[T99l. [1991 ISteidel et al.l 
119951 [2000f ). The color selection applied in searching for 
LBGs requires bright rest-frame UV continua and blue 
colors, therefore selecting galaxies with recent star forma- 
tion and stellar populations dominated by young stars. 
LBGs have greatly impacted investigations of galaxy evo- 
lution. For example, significant effort has been devoted 
to studying the cosmic star formation history and partic- 
ularly in measuring the peak (and decline at z > 3 — 4) of 
the gl obal star form a tion rate (SFR) densi ty using LBGs 
(,Madau et all [l99l IBouwens et all [2001 [2009 . .20111 
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However, dust attenuation in UV-selected galaxies poses 
challenges to measuring the total SFR in these galaxies. 
Other samples o f galaxies {e.g., submillimeter selected; 
iBlain et al.lll999l ) have been studied to provide additional 
constraints on the SFR density. 

X-rays offer another way to explore the dust- 
unobscured cosmic star formation history. X-rays in 
"normal galaxies" {i.e., not dominated by active galac- 
tic nuclei, AGN) mainly originate from accreting X-ray 
binaries (XRBs, including ultraluminous X-ray sources). 
The 2-10 keV emis sion in normal local galaxies scales 
with S FR ( Ranalli eTall I200I iPersic k Rephaelil 120071 : 
iMineo et al.. 2012). However a number of other sources 
can contribute to the X-ray emission in normal galax- 
ies: supernovae and their remnants, and hot gas from 
starburst-drive n winds and outflows (see, e.g., review by 
lFabbianolll989f l. Low-luminosity AGN activity may also 
contribute to the X-ray luminosity in normal galaxies. 

Within the XRB population, high mass X-ray binaries 
(HMXBs) are short-lived, tracing recent star formation 
activity (on timescales ~ 10^^^ yrs), while low mass X- 
ray binaries (LMXBs) trace older stellar populations (for 
timescales > 10*~^ yrs) and the 2-10 keV X-ray lumi- 
nosity scales with the stellar mass (M*) of the galaxy. 
Assuming the following analytic parameterization for lo- 
cal (within a distance of 60 Mpc) normal galaxies: 

ix,2-iokcV - ix(LMXB) + Lx(HMXB) (1) 
= aM, +/?SFR (2) 

iLehmer et al.l ()2010D measured constants of a 



(9.05 ± 

0.37) X 10^8 erg s'^ M'^ and /3 = (1.62±0.22) x 1 0=^^ erg 
s-i (M0 yr-i )-i (see also lColbert"eraL|[25ol . As a 
relatively dust- insensitive probe of past and present SFR, 
the average X-ray properties of star- forming galaxies give 
an independent measurement of the cosmic star forma- 
tion history. In addition, the wide range of evolutionary 
timescales, SFRs, and metallicities probed by deep sur- 
veys to high redshift offer constraints on binary evolution 
th eorie s {e.g., Fragos ct al. 2012, hereafter, referred to 
asEH). 

Several studies have shown the utility of stacking 
deep X-ray data at the optically determined locations 
of star-forming galaxies to measure the average X- 
ray properties of distant LBG populations that are 
individually undetected. Averaging the X-ray counts 
at the known positions of the targets enhances the 
signal-to-noise for source populations with fluxes be- 
low the detection limit. Through stacking, the re- 
lationship between X-rays and other galaxy proper- 
ties {e.g., SFR, stellar mass, dust attenuation) have 
been inves tigate d for low redshift galaxie s (z < 1.4; 



Ptak ct al. 2001; Hornschemeier et ar."2002'; 'Lair d et al.l 
2005; Lehmer et al. .2007, ,2008, : .Watson et al., 1200 
Svmeonidis et aLl " 201 ll). intermediate redshift LBGs 
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200. 



< z < 3 
and distant 



Brandt et all [200l iReddv fc Steidei 
LBGs (3 < z < 6; iBrandt et al 
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2001': 'Laird et al. 2006l: rLehmer et al.ll2005bHCowie et al 
2011). These stacking studies highlight the use of X-rays 
from normal galaxies to study t he star formation history 
of the Universe. For example, iSvmeonidis et al.l (|2011[ ) 
correlate X-ray emission with infrared (IR) emission for 
galaxies at < z < 2 and show that X-ray emission 
can be calibrated as a useful star formation indicator to 
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Fig. 1. — We show the areal coverage of the different LBG sam- 
ples, overlaid on the 4 Ms Chandra full band (0.5-8 keV) image 
of CDF-S. The F225W- (z = 1.5), F275W- {z = 1.9) and F336W- 
{z =2.5) dropouts fall within the blue region; U-dropouts (2 = 3) 
lie within the green region; B435- (z = 3.8), Vgoe- = 5.0), and 
1775- {z = 5.9) dropouts are located within the orange contour; 
Z850- (z = 6.9) and Y105- (2 = 8.0) dropouts are bounded by the 
red dashed lines. Additional Vgoe-j i775-> zgso-, and Yio5-dropouts 
reside in the CANDELS fields (shown in gray and white dashed 
lines). The thick black line shows the Chandra area with off-axis 
angle <7'. 



measure the global SFR density. 

In this paper, we perform stacking analyses of LBGs, 
using the deepest X-ray obser vation s to date (4 Ms Chan- 
dra observations of CDF-S; iXllf ). adding newly dis- 
covered LBGs and covering the broad redshift range, 
1-5 ^ z < 8 {i.e., correspondi ng to when the age of the 
Universe was 4.5-0.6 Gyrs). ICowie et al.l (|201l[ ) have 
completed a study of X-ray emission in this sample to 
find that X-ray emission from z = — 8 LBGs is con- 
sistent with normal galaxy populations, provided that a 
reasonable range o f dust attentiation factors (~ 3 — 5) are 
present. However. ICowie et al.l (|2011[ ) did not attempt to 
correct for dust attenuation and search for low levels of 
evolution in the X-ray emission. For the first time, we 
measure an empirical law that relates X-ray luminosity, 
SFR and redshift (z) and con ipare our results with X-ray 
binary evolution models feee lF13 ). 

Our paper is organized as follows: we introduce the 
samples and describe our stacking analysis techniques 
in Section [21 we discuss the individually detected X-ray 
sources in Section 13.11 and the stacking results for z < 4 
and z > 5 in Sections 13. 21 - [3^ finally, we offer our con- 
clusions in Section 21 The Galact ic column densitie s are 
8.8 X 10^9 cm-2 for the E-CDF-S (|Stark et al.lll992n . All 
of the X-ray fluxes and luminosities quoted throughout 
this paper have been corrected for Galactic absorption. 
We assume the standard ACDM cosmology: r2A/=0.3, 
r2A=0.7, and Ho=70 km s~^ Mpc~^, and th e initial mass 
function (IMF) described bv lKroupal (|200H) . 

2. DATA AND ANALYSIS 
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TABLE 1 

Lyman Break Galaxy Samples for Deep Chandra X-ray' Study 









Area 




(muv) 


(SFRuV,corr> 


reference 


z 


dropout band 


instrument 


(arcmin^) 


NlBGs 


(mag) 


(Mo yr-i ) 





1.5±0.4 


F225W (2372A) 


HST /WFC3 


50 


48 


24.5 


13 


Oesch et aL C2010a)t 


1.9±0.4 


F275W (2710A) 


HST /WFC3 


50 


91 


24.8 


30 


£csch ct aL (2010a) t 


2.5±0.6 


F336W (3355A) 


HST /WFC3 


50 


359 


25.4 


17 


Oesch et aL ('2010a)t 


3.0±0.2 


u (3570A) 


CTIO + HST /ACS 


160 


361 


25.3 


33 


Lee et aL C2006~|t 


3.8±0.3 


B435 (4317A) 


HST /ACS 


196 


2098 


26.5 


8 


Bouwens et aL (2007) 


5.0±0.3 


V6O6 (5918A) 


HST /ACS 


196 


445 


26.5 


10 


Bouwens et aL (20071 




i775 (7693A) 


HST /WFC3 (CANDELS) 


105 


768 (700)* 


26.8 


10 


Bouwens et aL (2012') 


5.9±0.3 


HST /ACS 


196 


181 


27.0 


5 


Bouwens et aL (2007') 




Z850 (9055A) 


HST /WFC3 (CANDELS) 


105 


218 (208)* 


27.0 


8 


Bouwens et aL (2012') 


6.8±0.4 


HST /ACS 


150 


73 


27.9 


4 


Bouwens et aL (2011) 






HST /WFC3 (CANDELS) 


105 


41 (31)* 


26.7 


11 


Bouwens et aL (2012') 


8.0±0.5 


Y105 (l.OSS/xm) 


HST /WFC3 


150 


60 


28.0 


2 


Bouwens et aL (2011) 






HST /WFC3 (CANDELS) 


105 


24 (22)* 


27.0 


5 


Bouwens et aL (2012) 



^ Mean apparent rest-frame UV magnitude. 

iVIean dust-corrected UV SFR. 
^ catalog of sources was acquired via private communication. 

The total number of LBGs in the dropout sample, with the number of unique (not included in ACS sample) LBGs in parentheses. 



For our analysis of the X-ray properties of LBGs, we 
use the pubhcly available 4 Ms CDF-S catalog and data 
products from XI f^. 

We use 1-5 < z < 8 LBG samples compiled from 
sever al catalogs ( see Ta ble [T]): z ~ 1.5, 1.9 and 2 5 
from lOesch etall (|2010al) . z 3 fromlLee et all ([200l . 
z ss 3.8, 5.0 and 5.9 from Bouwen s et al and 
z K. 6.8 and 8.0 from Bouwens ct all ()2011h . We in- 
clude additional z 5 — 8 LBGs from the Multi-Cycle 
Treasury Program Cosmic Assembly Near-Infrared Deep 
Extragalactic Legacy Survey (CANDEL S) CDF-S field 
(jBouwens et al.ll2012l: lOesch et al.ll2012f ). We refer read- 
ers to these papers for details about the LBG selection, 
however the basic technique for selecting these galax- 
ies involves using photometric color selection to isolate 
galaxies with relatively flat spectra that drop out of the 
bluest filter, while avoiding the expected color tracks of 
nearby interlopers (stars, low-redshift elliptical galaxies). 

The areal coverages of the LBG samples are shown 
in Figure [Don top of the Chandra 0.5-8 keV full-band 
image. The blue region marks the LBGs selected from 
the Hubble Space Telescope {HST) Wide Field Cam- 
era 3 (WFC3) Early Release Science data, i.e., galax- 
ies dropping out of the F225W (z « 1.5; 48 galaxies), 
F275W(z « 1.9; 91 galaxies), and F336W(z « 2.5; 
359 galaxies) filters. The green region contains the 361 
z « 3.0 LBG sample, selected as U-dropouts, based 
on ground-based U-band observations taken on the 4m 
telescope at Cerro Tololo Inter-American Observatory 
(CTIO) of the Great Observatories Origins Deep Survey- 
South (GOODS-S) field. The orange contour bounds the 
samples of 2098 B435-(^ ~ 3.8), 445 Veoe- {z « 5.0) 
and 181 1775- (z sa 5.9) dropouts selected using the Ad- 
vanced Camera for Surveys (ACS) instrument aboard 
HST, which includes the GOODS-S field and two ad- 
ditional regions from parallel fields for NICMOS Ultra- 
deep Field (UDF; PI: Roger Thompson: iThompson et all 



Data products are available online at 

|http : //www2 ■ astro . psu . edu/- {}niel/cdf s/cdf s- Chandr a . html | 



120051) . The red dashed regions mark the locations of 73 
z « 7 and 60 z w 8 LBGs, selected from dropping out 
of the Z850- and Y105 filters using HST ACS and WFC3 
infrared data. The CANDELS fields, shown with gray 
and white dashed lines (CANDELS-Deep in upper re- 
gion; CANDELS- Wide field in lower), add 700 Vgoe-, 208 
1775-5 31 zgsQ- and 22 Yio5-dropouts. This information is 
summarized in Table [T] 

Contamination from various sources {e.g., low-mass 
stars, AGN, low redshift galaxies, transient sources, pho- 
tometric scatter and spurious sources) in the LBG sample 
is expected to range from 3% t o 10% , an d is discussed in 
more detail in iBouwens et al.l (|2007[ ) and iBouwens et al.l 

ipoTl 

2.1. LBG samples: Determining UV-based 
Dust- corrected SFRs 

In order to study the relationship between X-ray lu- 
minosity and SFRs in LBGs, we estimate their SFRs 
based on rest-frame U V properties. We use the relation 
given m Equation 1 of iBell et al.l (pOOl to convert rest- 
frame UV luminosities into UV-derived SFRs, given by 
the right-side of this equation: 

SFR (Moyr-i) = 9.8 x 10-"(Lir + 2.2Luv)/L0 (3) 

where Lyv is the UV luminosity (= i'li/,2800i where v is 
the frequency, and lu.2s,Q0 is the monochromatic luminos- 
ity measured at rest-frame 2800A) and L© is solar lumi- 
nosity (L0 = 3.84 X 10"^'^ erg s~^). Based on the available 
LBG catalogs, we note that the given UV magnitudes 
correspond to rest-frame 1500A for some galaxy samples 
(z sal. 5, 1.9, and 2.5) and to 2800A for the other sam- 
ples. Since the LBG galaxy spectrum (L^) is expected 
to be very flat in this spectral region there are negli- 
ible corrections to the SFR based on these differences 
Kennicuttlll998D . 

We derive dust corrected SFR s follow ing the pre- 
scription given in IBouwens et al.l (|2011[ ). Using the 
measured UV-continuum slopes (/3) in z « 2 — 
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TABLE 2 

Stellar Mass Properties for 2 < 6 Lyman Break Galaxies 



(logM,) Med (log M.)* (logSFR) Med (log SFR)* (logsSFR) Med (log sSFR)* 
N^.t.h^' Ntot'' (Mq ) (Mq ) (M0 yr-i ) (Mq yr"! ) (yr^i) (yr-i) 



All SFRs 



1.5 


42 


47 


9.6 ±0.6 


9.5 


1.0 ±0.3 


1.0 


-8.6 ±0.6 


-8.4 


1.9 


59 


91 


9.5 ±0.5 


9.4 


1.4 ±0.3 


1.4 


-8.1 ±0.5 


-8.1 


2.5 


127 


358 


9.4 ±0.7 


9.4 


1.2 ±0.4 


1.2 


-8.2 ±0.6 


-8.2 


3.0 


156 


360 


9.5 ±0.6 


9.5 


1.2 ±0.3 


1.1 


-8.3 ±0.5 


-8.3 


3.8 


982 


1847 


8.9 ±0.6 


8.9 


0.8 ±0.4 


0.7 


-8.1 ±0.5 


-8.1 


5.0 


414 


874 


8.9 ±0.8 


8.9 


0.8 ±0.3 


0.8 


-8.1 ±0.7 


-8.1 


5.9 


136 


292 


8.9 ±0.9 


9.0 


0.7 ±0.3 


0.7 


-8.2 ±0.8 


-8.3 


6.8 


4 


65 


9.5 ±0.6 


9.9 


1.5 ±0.2 


1.5 


-8.0 ±0.4 


-8.0 



medium SFRs : SFR= 5 - 15 Mq yr^^ 



1.5 28 31 9.6 ±0.6 9.6 1.0 ±0.1 1.0 -8.7 ±0.7 -8.5 

1.9 17 28 9.6 ±0.8 9.3 1.1 ±0.1 1.1 -8.5 ±0.8 -8.3 

2.5 49 167 9.2 ±0.7 9.2 0.9 ±0.1 0.9 -8.3 ±0.7 -8.2 

3.0 89 224 9.3 ±0.5 9.3 1.0 ±0.1 1.0 -8.3 ±0.5 -8.3 

3.8 357 580 9.0 ±0.5 9.1 0.9 ±0.1 0.9 -8.1 ±0.5 -8.1 
5.0 196 348 9.0 ±0.8 9.0 0.9 ±0.1 0.9 -8.1 ±0.7 -8.1 

5.9 53 80 8.9 ±0.9 9.0 0.8 ±0.1 0.8 -8.0 ±0.8 -8.0 



high SFRs: SFRs = 15 - IOOMq yr^^ 



1.5 8 9 9.8 ±0.5 9.8 1.5 ±0.1 1.5 -8.3 ±0.6 -8.1 

1.9 40 61 9.5 ±0.4 9.5 1.5 ±0.2 1.5 -8.0 ±0.3 -8.0 

2.5 56 105 9.7 ±0.5 9.7 1.5 ±0.2 1.5 -8.2 ±0.5 -8.1 

3.0 64 119 9.8 ±0.5 9.7 1.5 ±0.2 1.4 -8.3 ±0.5 -8.2 

3.8 156 231 9.5 ±0.5 9.5 1.4 ±0.2 1.4 -8.1 ±0.5 -8.1 
5.0 57 112 9.5 ±0.7 9.5 1.4 ±0.2 1.3 -8.1 ±0.7 -8.2 

5.9 10 14 9.9 ±0.4 10.0 1.3 ±0.2 1.4 -8.6 ±0.3 -8.5 
6.8 4 9 9.5 ±0.6 9.9 1.5 ±0.2 1.5 -8.0 ±0.4 -8.0 



Median values 

^ The number of matches with lXue et al.l I I2012I) with \Az\ < 0.5 
^ The total number of LBGs within lXue et al.l 120121) footprint. 

7 LBGs f Bouwensejal] [2009t IBouwens et all [20TTI ). 
IBouwens"et al. (2 0111 ) provide a conversion of the ob- 
served UV absolute magnitude and redshift into an es- 
timated p. We apply the relation b e tween (3 and FUV 
dust attenuation from iMeurer et all (|1999[ ) to estimate 
the extinction factor for each galaxy, which also takes 
into account the typical scatter of /3 values 0.36) at 
an observed UV magnitude and redshift. These extinc- 
tion factors range from 3-6 for the z < 3 LBGs and are 
<2 at z > 3, with typical uncertainties of ^ 0.3 (see 
Table |4]) . As an additional check on the dust-corrected 
SFRs, we compare the average dust-corrected UV SFR 
with radio- and UV+IR-derived SFRs for each redshift 
sample in Section [2.51 

We show the SFR distributions of LBG samples in Fig- 
ure O We divide our LBG samples into separate SFR 
bins (marked by the dotted lines): SFR/(M0 yr^^)= 
1-5 (low SFR), SFR/(M0 yr-i)= 5-15 (medium SFR) 
and SFR/(M0 yr-i)= 15-100 (high SFR) with the num- 
bers of galaxies that fall in each bin labeled. This choice 
of SFR binning was selected to match th e SFR binning 
at low redshifts (|L08l : iLehmer et al.|[20Tol) for direct com- 
parison. The dark gray arrows mark the mean for each 
SFR bin and the light gray marks the expect mean, cal- 
culated using the UV luminosity function (UVLF; the 
form is described by the parameters given in the papers 
listed in Table [T] for each LBG sample) to predict the 
number of galaxies beyond the observed limit. The z ^ 3 
medium SFR samples appear to suffer from incomplete- 



ness, in that the observed mean SFRs are overestimated 
compared to the true mean for LBGs. We discuss the im- 
plications of this incompleteness on our results in Section 
[3J and 1331 

To compar e stellar m a sses, M», for these LBGs, we 
match to the iXue et all ()2012[ ) catalog, which uses the 
available (12 bands) photometric data to estimate stellar 
masses. We only consider ma tches where the redshifts 
(between our sample and the iXue et all 120121 catalog) 
agree within |Az| < 0.5. The number of matches de- 
creases for fainter and higher redshift objects. Only 4 of 
the z > 6 LBGs were mat ched to sources in the lXue et al.l 
(|2012l ) catalog; since the IXue et al.l (|2012[ ) catalog has a 
S-c Z85o-band limiting magnitude of 28.1, it is likely that 
most of ou r highest redshift LBGs are undetected in this 
filter. The lXue et al.l ()2012[ ) catalog gives the most likely 
redshift, but it is possible that the second most likely red- 
shift would match the LBG selection redshift {e.g., see 
Section O about XID=371). While the |Az| < 0.5 cri- 
terion diminishes the sample size significantly, the av- 
erage and median offer valuable information about 
the samples: (logM*) — 8.9 — 9.6 for the full sample; 
8.9—9.6 for the medium SFR sample; and 9.5—9.8 for the 
high SFR sample. The specific SFRs (sSFRseeSFR/M,) 
are (logsSFR) = —8.7 to —8.0. The median and mean 
values for M*, SFR, and sSFR are given for each sub- 
sample (redshift and SFR range) in Table |2| Given the 
high sSFRs for these galaxies, the stellar populations in 
these galaxies are likely relatively young, exhibiting high 
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Fig. 2. — We show the dust corrected UV-derived SFR distribu- 
tions of the LBG samples used in this study. The dotted lines 
mark our SFR cuts for two separate stacking bins: SFR/(Mq 
yr-l)= 5-15 and SFR/(M0 yr-i)= 15 - 100 with the num- 
bers of galaxies that fall in each bin labeled. The arrows mark the 
observed mean (dark gray) for the different SFR bins compared 
to the expected mean (light gray; based on the UVLF). The low- 
est SFR bins do not yield any stacked detections and suffer from 
incompleteness; therefore, they are not discussed further. 

present to past SFRs. 

2.2. LBG samples: individually detected X-ray sources 

Since our primary goal is to study the X-ray emission 
arising from star formation in these high-redshift LBGs, 
we avoid including AGN in our X-ray stacking analyses. 
We assume that sources that have been individually de- 
tected at these redshifts (and therefore with Lx > 10^^ 
e rg s~ ^) must be dominated by emission from AGN. 

Kill have matched the X-ray detected sources to 
optical, infrared and radio catalogs using probabilis- 
tic matching, which estimates the false match statis- 
tics by applying a Monte Carlo technique (see also 
iBroos et al.l 120111) and accounts for positional uncer- 
tainties (whi ch va ry across the Chandra field with off- 
axis angle). IX 111 found an offset between the the op- 



tical/infrared positions and radio positions (the X-ray 
positions are fixed to the radio frame): RAx. radio = 
RAopt.iR+0.175"and Decx,radio = DeCopt,iR-0.284". 
The coordinates (shifted to the radio/X-ray system) for 
the best matched GOODS-S or GEMS sources are pro- 
vided, which we match to our LBG positions within 0.5". 
We find 20 X-ray detected LBGs (which are likely AGN) 
and tabulate their properties in Table [S] 

In this Table, the redshifts based on the LBG sample 
selection (z) are listed along with the redshifts from the 
1X1 II catalog, based on spectroscopic (s) or other photo- 
metric data (p). The effective photon index (F) was de- 
termined for all sources with >200 full band ( 0.5"8 keV) 
counts using basic spectral fitting (described in IX 111 ) and 
assuming an absorbed power law. The last column in this 
table includes notes about previ ously detected and stud- 
ied so urces from CDF-S 1 M s KM) . E-CDF-S catalog 
(L05a: [Lehmer et"al] I2005a[) . and radio (R) detections 
based on our VLA data analysis (see Section 1^31. T hese 
sources are discussed in more detail in Section [5T1 



2.3. X-ray stacking of individually undetected LBGs 

While even 4 Ms of exposure time is not sufficient for 
detecting X-rays from individual star-forming galaxies 
at z > 1 (with the possible ex c eption of submilli meter 
galaxies; lAlexander et al.l 120051: iLaird et all I2O1O0 . the 
X-ray emission from large numbers of galaxies can be 
stacked to provide the average X-ray properties of the 
stacked sources. We fo l low t he stacking technique out- 
lined in ILehmer et al.l ()2008l hereafter, LOS) to stack 
LBGs in bins of SFR (shown in Figure [5]) and z. 

First, we select so urces that have no X-ray detections 
in the catalog from IXIII within twice the 99% encir- 
cled energy fraction (EEF) radius (note that this value 
changes with off-axis angle) in order to avoid contamina- 
tion from the Chandra point spread function (PSF) from 
detected sources. We only include sources with off- axis 
angles < 7' from the central pointing (thick black circle 
in Figure [T]) because the increasing PSF at large off-axis 
angles degrades the X-ray sensitivity significantly. These 
two criteria eliminate 856 and 837 LBGs, respectively, 
from our original sample of LBGs. 

Rec ent stacking results of z > 6 LBGs bv lTreister et al.l 
(|2010| ) have generated discussion about proper back- 
ground subtraction in order to avoid artificial stacked 
detections {e.g., by setting the clipping threshold of the 
background pixels too low : seCiCowie et al. 2011; Willot3 
1201 If) . We refer readers to ILehmer et al.l (|200"5a) and LOS 
for detailed descriptions about our background subtrac- 
tion technique using background maps, which do not in- 
clude any clipping of photons. Briefly, these background 
maps were created by filling the regions near detected 
sources (twice the 90% PSF EEF radius, rgo) with noise, 
which is generated based on the probability distribution 
of counts in the local background (within an annulus with 
radii between 2 and 4 x rgo). 

Throughout the remainder of this paper. X-ray lumi- 
nosities refer to the rest-frame hard (2-10 keV) energy 
band. These luminosities are k-corrected from the ob- 
served soft band (0.5-2 keV) since this band is closer 
to the rest-frame 2-10 keV band, and therefore the k- 
correction is smaller and less affected by the assumed 
spectral shape {i.e., photon index, F). For example, at 
the median redshift of our LBG samples (z sa 3) the 
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TABLE 3 

Individually Detected X-ray Sources Associated with Ly'man Break Galaxies 



log Luv"" SFRi^ O.S-SkeV^ log Lo.5-8'' 



XID 


RA 


Dec 


2 


^Xuo 


(erg s ^) 


(Mq yr 1 ) 


(counts) 


(ergs ^) 


r 


Notes'' 


344 


53.10485 


-27.70521 


1.5 


1.6 


s 


43.91 


6.3 


1811.3150-? 


43.79 ±0.01 


2.0 ±0.1 


L05a, 


R 


405 


53.12284 


-27.72279 


1.9 


1.6 


V 


43.97 


13.4 


59.0l}5.i 


42.93 ±0.11 


< 0.1 


R 




308 


53.09392 


-27.76774 


1.9 


1.7 




43.95 


12.3 


674.5l^j;;g 


43.92 ±0.02 


0.2 ±0.1 


L05a 




326 


53.10082 


-27.71601 


2.5 


2.3 




44.58 


54.9 


234.4+?nI 
— 20.5 


43.46 ±0.04 


1.6 ±0.2 


R 




137 


53.03334 


-27.78257 


3.0 


2.6 




44.36 


18.9 


1249.4+tS-o 


44.45 ±0.01 


1.3 ±0.1 


R 




490 


53.14880 


-27.82112 


3.0 


2.6 




44.31 


16.4 


84 6+1^ 

" ^ — 10.9 


43.43 ±0.06 


0.4 ±0.2 


R 




386* 


53.11792 


-27.73432 


3.0 


3.6 




44.38 


20.1 


37 4+"'! 
' 10.6 


42.99 ±0.13 


> 1.0 


R 




254 


53.07600 


-27.87816 


3.0 


2.8 




44.67 


44.1 


246 2+?g1 


43.65 ±0.04 


1.8 ±0.2 


L05a, 


L05 


563 


53.17439 


-27.86735 


3.0 


O.D 


s 


45.43 


365.0 


1389 l+lf5 

-Ll_J(_Jt/. -L 41.4 


44.51 ±0.01 


1.4 ±0.1 


ijuoa 




573* 


53.17848 


-27.78403 


3.0 


3.2 


s 


44.31 


16.4 


685.1129.? 


44.14 ±0.02 


1.4 ±0.1 


L05a, 


L05 


577 


53.18015 


-27.82060 


3.0 


1.9 


s 


45.41 


342.6 


2437.0l53;3 


44.67 ±0.01 


1.6 ±0.1 


L05a, 


L05 


588 


53.18464 


-27.88092 


3.0 


3.5 


s 


44.80 


64.0 


322.4l24;3 


43.79 ±0.03 


1.8 ±0.2 


L05a, 


L05 


388* 


53.11858 


-27.88478 


3.8 


3.0 


p 


44.12 


7.2 


41 3+12.7 


43.38 ±0.13 


< 0.3 






374 


53.11163 


-27.86078 


3.8 


3.7 


V 


44.48 


19.3 


< 26.5 


< 43.11 


< 0.7 






262 


53.07848 


-27.85984 


3.8 


3.7 


s 


44.50 


20.2 


214.6ll?:«, 


44.19 ±0.04 


0.1 ±0.1 


R 




150* 


53.03989 


-27.79846 


3.8 


3.3 


p 


44.11 


6.9 


38.6llJ;« 


43.32 ±0.13 


< 0.5 






100* 


53.01660 


-27.74484 


3.8 


3.9 


p 


44.42 


16.1 


93.4ll^;3 


43.63 ±0.09 


> 1.1 






546 


53.16528 


-27.81405 


3.8 


3.1 


s 


44.75 


40.0 


1231. Olgg;? 


44.84 ±0.01 


0.4 ±0.04 


L05a, 


L05, 


458 


53.13850 


-27.82112 


3.8 


3.6 


p 


44.49 


19.8 


27.6l^-? 


43.26 ±0.12 


0.2 ±0.04 






371 


53.11158 


-27.76777 


5.0 


3.1 


p 


44.16 


7.0 


173.311^;;^ 


44.19 ±0.04 


0.6 1°;? 







^ Observed {i.e., not dust-correeted) UV luminosity. 

^ Dust-corrected UV-derived SFR. 

^ Counts refer to the aperture-corrected (background-subtracted) net counts for the full-band (0.5 — 8 kcV) from [xTU 

*^ Lo,5 — 8 is the full band luminosity derived using the 0.5—8 keV observed flux, k-corrected by (1 -|- 2:)^^^, assuming V — 1.8 (typical photon 
index for AGN), and z (not ^xuc)- Upper limits correspond to 2.5-(T limits for sources that were detected in only the hard (2—8 keV) band, while 
undetected in the soft (0.5 — 2 kcv) and full bands. 

*^ The code L05 refer to X-ray detected sources in and were consequently left out of their stacking analyses of Chandra 1 Ms data. L05a 

refers to the X-ray sources identified in tire E-CDF-S catalog JLehmer et al.ll2005al ). R indicates that this source was detected in the VLA radio 
data with S/N>3 (see Section [231 . 

Found in a close pair (|A2:| < 0.5 and angular separations of 3— 4" (corresponding to physical separations of 20—30 kpc at 2 ^ 3 — 4) with another 
LBG. Sec discussion in text fSection l3.lt . 



observed 0.5-2 keV band corresponds to rest-frame 2-8 
keV. The fc-correction can be calculated as follows: 
p(2-r) _ p(2-r) 

, _ ^out,2 ^out,l ~Vr-2.0) 

l^corr - (2-r) _ p(2-r) ^ > 

= 2.472 + (5) 

using r = 1. 7, based on the ex pected photon index for X- 
ray binaries (jPtak et al.lll999f ) and in the input and out- 
put energies: Eout,i =2 keV, Eout,2 =10 keV, E^n,i =0.5 

keV, =2 keVH 

We consider two sources of error in our stacked mea- 
surements. The first, CTp, is Poisson noise and calculated 

simply as a/B, where B represents background counts. 
This is the error that we use to determine the signal 
to noise. In addition, we also calculate the errors from 
boot-strapping (crboot)i which measure how the contribu- 
tion of individual sources may affect the average value. 
To determine our boot-strapping errors, we randomly re- 

1^ We note that we attempted to stack in the observed 2-8 keV 
bandpass, however we did not obtain > 2(t detections in any cases. 
Ideally, we could use the 2-8 keV/0.5— 2 keV stacked band ratios 
as proxies for the spectral slopes of our stacks. Unfortunately, the 
upper limits of the band ratios are not well constrained (Fefi >0.3- 
0.7) for all cases. Therefore, we limit further discussion to results 
based on stacking the observed soft 0.5-2 keV band. 



sample the LBGs for each stacking bin a large number 
of times (5000) and repeat our stacking analysis on these 
samples. The random resampling will duplicate some 
values while eliminating others each iteration to statisti- 
cally quantify the effect of individual sources on the stack 
(~ 37% of the v alues are replaced by duplicated values; 
see lEfronlll982l ). The boot-strapped errors refer to the 
standard deviation of the stacked values from this ran- 
dom resampling. Bootstrap errors include the variation 
due to background fluctuations (Poisson error), therefore 
the uncertainties quoted in our analysis are the bootstrap 
errors. 

To study the relationship between X-ray luminosity 
and SFR of z < 4 LBGs, we stack the LBG samples in 
two SFR bins: SFR/Mq yr"! = 5-15 and SFR/M© 
yr"^ = 15—100. For the remainder of our study, we 
require a S/N cut of 2.5-cr for a stacked detection (which 
corresponds to >99.38% confidence limit for one-sided 
gaussian distributions). Since the low SFR samples do 
not yield any detections and suffer greatly from incom- 
pleteness, we do not include them in our discussion. The 
results of our stacking analyses are discussed in detail in 
Section O 



2.4. Contribution of AGN to average X-ray properties 
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TABLE 4 

Comparing SFRs: UV, Radio, and IR 



41 42 43 

log 2-10 keV Luminosity 



44 



Fig. 3. — We show the cumulative AGN fraction vs. 2-10 keV X- 
ray luminosity for SFR=5— 15 Mq yr~^ (medium SFR; open green 
symbols) and SFR=15-100 Mq yr"! LBGs (high SFR; filled blue 
symbols) and for z < 1.4 star-forming galaxies from LOS (gray 
triangles; filled for high SFR and open for medium SFR). The 
different symbols describe different LEG samples, as listed in the 
legend. The solid lines mark the error-weighted exponential best 
fit to the combined dataset for the high SFR (blue) and me dium 
SFR (green) data; best-fits determined from Equation 8 in IL081 
using mean SFRs from the full sample for high SFR and medium 
SFR galaxies arc shown as blue dashed and green dotted lines, 
respectively. 

The 20 X-ray detected sources in Table [31 which were 
classified as AGN following the criteria outlined in Sec- 
tion l2.21 are useful for assessing the level of AGN contam- 
ination in the stacked samples discussed in the previous 
section f§ 12.3 1. 

In Figure [21 we show the cumulative AGN fractions as 
a function of X-ray luminosity for high SFR (15 — 100 Mq 
yr~^ ; filled blue symbols) and medium SFR (5 — 15 Mq 
yr~^ ; open green symbols) LBG samples. The different 
redshift samples are displayed with different symbols, as 
described in the legend. 

The data are not sufficient to determine whether there 
is any variation in the AGN fraction with redshift. Fo- 
cusing on the high SFR LBG sample (filled blue sym- 
bols), where we have some redshift information, there 
is no obvious redshift trend {e.g., the z < 3 point ap- 
pears consistent with the z « 3.8 data, while the z sa 3 
values are higher). We proceed by including high SFR 
data from z < 1.4 star-formi ng ga laxies from iLOSl (filled 
gray triangles) . The fit from 'L08' (blue dashed line) , us- 
ing the mean SFR from the high SFR sample, compares 
well with our error-weighted fit (solid blue line) to all 
the dat a (iLOSi and our high SFR sample points). The 
fit from lLOST Tgreen dotted line; open gray triangles show 
medium SFR galaxies from this study), using the mean 
SFR from the medium SFR sample, appears consistent 
with our medium SFR fit (dotted green line) to data 
(open green circles). 

Our higher redshift data probe a larger volume and 
extend the relatio n to higher Lx- At lower Lx, the 
difference between ILOSl may be from incompleteness at 
these lower X-ray luminosities. Therefor e we estimate 
the AGN contribution from Equation 8 of IL08I and sub- 
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61 
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Q +6 
-2 


< 713 


<2200 


8.0 


44 


r, +3 

^ -1 


+3 
^ -1 


< 1064 
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The mean UV, dust-corrcctod UV, radio and UV+IR SFRs for full 
sample of LBGs (only excluding individually-detected X-ray sources 
and those sources within twice the 99% EEF radius from other X-ray 
sources) . 

tract this from the stacked X-ray luminosity as described 
in LOSl Since the AGN fractions from our analysis appear 
lower (or equal) to those from lL08l we err on the side 
of over-compensating for the AGN contribution. Nev- 
ertheless, these corrections are smaller than the errors 
on Lx (typically ~20%). We note that ILOSl estimated 
50 — 70% AGN contamination in z w 3 LBGs, based 
on CDF-S 1 Ms analysis. Our increased depth, which im- 
proves identification and removal of AGN, may account 
for the significant decrease in our measured AGN frac- 
tion; since low luminosity AGN are more numerous, sur- 
vey depth does not necessarily scale proportionally with 
AGN fraction. 



2.5. Comparing average SFRs using multiwavelength 
data: radio and IR SFRs 

We perform stacking of radio and IR observations to of- 
fer independent tests for estimating unattenuated SFRs 
for th e z < 5 LBG s amples. Us ing VLA ( Miller et al.l 
[2(M _in prep; see also [Miller et al.. ,2008') and FAGS Her- 
shel (jElbaz et alll201lD data for the CDF-S, we have 
stacked all of the LBGs that are not X-ray detected. 
These measurements provide consistency checks between 
the radio-derived, UV-I-IR SFR, and the dust corrected 
UV-derived SFRs. 

We generated 50" x 50" cutouts of the radio image 
around each LBG position, and proceeded to stack these 
cutouts for each LBG redshift sample. These LBG sam- 
ples exclude the same galaxies as the X-ray stacking sam- 
ple {i.e., those sources that were individually X-ray de- 
tected or near, within twice the 99% EEF radius, other 
X-ray detected sources). Incidentally, 8 of the 20 in- 
dividually X-ray detected LBGs (listed in Table [3]) were 
also individually detected (with > 3a) in the VLA obser- 
vations (marked "R" in the last column). Two stacking 
procedures were used: the first used a straight median for 
each pixel in the stack, and the second used an average 
with rejection of the highest and lowest pixel. The lat- 
ter method is designed to prevent real sources within the 
50" fields from seeping into the average and hence out- 
put stack image, which is more likely when the number 
of images to stack is large because of the surface density 
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of such bright sources. These real sources are often ex- 
tended radio emission from unrelated sources, although 
a small number of the LBGs themselves were individu- 
ally detected (7 of the 13 radio-detected LBGs with radio 
S/N > 5(J were also X-ray detected and are marked with 
'R' in the last column of Table |3]) and the rejection pre- 
vented them from unduly influencing the resulting stack. 
For the z ~ 3.8 and z « 5.0 samples, the rejection was 
extended from the single highest pixel to the two highest 
pixels to ameliorate contributions of real detected sources 
and improve the cosmetics of the stacked image. 

The stacked images were then inspected to evaluate 
the resulting RMS noise level and any sources recovered 
by the stacking technique. The results of the stacks pro- 
duced by the median and average procedures were al- 
ways consistent in terms of the resulting RMS noise and 
for the sample that yielded a weak detection. The RMS 
of the stacks ranged from about 0.2/j,Jy to 1.3/iJy de- 
pending on the number of objects in the sample, down 
from the ~6.4^Jy noise level typical of the radio image 
across the CDF-S area. Stacks made from cutouts using 
arbitrary positional offsets demonstrated that the RMS 
noise values in the actual sample stacks were consistent, 
and these "blank sky" stacks produced no false detec- 
tions. The 1.4GHz radio luminosities were fc— corrected, 
assuming a standard power law with a spectral index 
of 0.7, and convert ed to SFRs using the relation given 
in lYun et al.l (j2001l see Equation 13; using the Salpeter 
IMF with mass limits of 0.1 and 100 Mq ), but divided 
by a factor of 1.7 to translate SFRs to the Kroupa IMF. 

Herschel-PAGS data covering the GOODS-S (see 
lElbaz et al.ll20ljl for a description of the observations) 
was downloaded from the HeDaM websiteE3. For the 
purposes of this paper we used only the 160/im data as, 
for our samples, it probes closer to the peak of the far- 
infrared SED than the other PACS wavebands and thus 
provides the most reliable indicator of the IR luminosity. 
To perform our stacking analyses we extracted 60pixel 
x 60pixel thumbnails (from both the science and RMS 
images) centered on the positions of all the X-ray unde- 
tected sources that are also undetected at 160/im. For 
each of our redshift bins we combined corresponding sci- 
ence thumbnails using a weighted mean, weighting by the 
\/y/{RMS) value at each pixel position. We use aper- 
ture photometry to measure the flux within a circular 
aperture placed at the centre of each stacked image. The 
error on the mean was calculated by randomly selecting 
(with replacement) two-thirds of each redshift sample, 
performing the our stacking procedure 1000 times, then 
calculating the standard deviation of the resulting flux 
distribution. Finally, to determine the significance of 
each stacked (mean) flux, we compared it against that 
obtained when stacking the same number of random po- 
sitions. The mean flux of the 160/im-undetected sources 
was then combined with the mean flux of the detected 
sources using a weighted mean (weighted according to 
the number of undetected/detected sour ces) to give the 
mean flux of all sources in that bin (see iMuUanev et al.l 
[20T1 for more details). We convert the stacked 160/im 
flux into a total IR ( 8-1000/im) luminosity using the tem- 
plate spectrum from lCharv fc Elbad ()200lD . Then, using 

^° [http : //hedam . oamp ■ f r/GOODS-Herschel/goods-south_data ■ php, 
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Fig. 4.— HST ACS zsso image shows XID=386, an X-ray de- 
tected LBG at z « 3, apparently interacting with another LBG 
(also at the same redshift). These LBGs have angular separations 
of 0.8", corresponding to a physical separation of 6 kpc. 




Fig. 5. — The observed 0.5-2 keV stacked images for the z a: 1.5 
and 4 LBGs with SFR= 15 - 100 Mq yr"! (left and middle) and 
z = 6 (right) LBG samples are shown. These images have been 
smoothed by a 3 pixel kernel using gaussian weighting. The z = 1.5 
image shows a 2.6-cr detection (near our detection threshold) and 
contains 6 galaxies (24 Ms exposure), z = A image shows a A.O-cr 
detection and contains 171 galaxies, corresponding to 684 Ms of 
exposure. The z = 6 image contains 220 galaxies (~ 880 Ms total 
exposure) with S/N= 1.0. The marked circles have 3"radii. 

Equation |3l where Lir is the total IR luminosity, we cal- 
culate the UV-I-IR SFRs. 

Table |4] summarizes the mean observed UV, dust- 
corrected UV, radio and UV-I-IR SFRs. Only the z = 2.5 
LBG sample was detected in the radio stacking at ^ 3a; 
none of the LBG samples were detected in the IR stacking 
analysis. Within errors, the dust-corrected UV-derived 
SFRs are consistent with the radio and UV-f IR SFRs. 

3. RESULTS AND DISCUSSION 

We divide our discussion of the results into three sec- 
tions. The first section includes a brief summary of inter- 
esting discoveries related to individually detected LBGs 
from the 4 Ms CDF-S catalog (jXllt ). The other sections 
relate to our stacking analyses. Based on our stacking 
results, we find that the z < A samples provide signif- 
icant stacked detections, whereas we were only able to 
measure upper limits for the z ^ 5 LBGs. Therefore, 
we consider it is most meaningful to split the stacking 
analysis into two parts: studying the evolution of the X- 
ray properties at z < 4 as related to X-ray binaries and 
star formation, and constraining supermassive black hole 
growth for z 3 — 8 LBGs. 
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TABLE 5 

Summary of Stacked X-ray/Star formation properties in z <4 LBGs 



z # of sources ttot" Net counts'^ S/N'^ log Lx'' log Luv log SFR<= log Lx/SFR 

(Ms) (ergs-i) (org s-i) (Mq yr^i) (erg s'l [M© yr-i]-i) 



5<SFR/(M0 yr-i )<15 



1.5 


21 


84 


40.3 ± 13.5 


3.6 


41.06 ± 0.13 


44.05 


1.01 


40.05 ± 0.13 


1.9 


15 


60 


< 23.3* 


2.0 


< 41.14 


43.94 


1.08 


< 40.06 


2.5 


99 


396 


< 60.3 


1.5 


< 41.06 


44.01 


0.95 


< 40.11 


3.0 


191 


764 


133.7 ± 35.5 


4.0 


41.20 ± 0.10 


44.11 


0.98 


40.21 ± 0.10 


3.8 


427 


1708 


< 123.6 


1.8 


< 41.07 


44.18 


0.93 


< 40.14 


15<SFR/(M0 yr-i )<100 


1.5 


6 


24 


15.3 ± 3.0 


2.6 


41.20 ± 0.08 


44.38 


1.46 


39.71 ± 0.08 


1.9 


28 


112 


48.9 ± 14.4 


3.8 


41.27 ± 0.11 


44.22 


1.48 


39.75 ± 0.11 


2.5 


62 


248 


81.1 ± 16.8 


4.3 


41.41 ± 0.08 


44.40 


1.51 


39.86 ± 0.08 


3.0 


104 


416 


90.0 ± 30.4 


3.7 


41.28 ± 0.13 


44.54 


1.50 


39.72 ± 0.13 


3.8 


171 


684 


126.2 ± 32.3 


4.0 


41.48 ± 0.10 


44.60 


1.44 


39.98 ± 0.10 



^ Total exposure time of the stack. 

^ Background-subtracted counts in observed 0.5-2 keV band. Errors include Poisson and bootstrap errors (as described in Section l2.3t . 

^ S/N is measured as S/V^, where S and B are net and background counts, respectively. 

^ Rest frame 2—10 keV luminosity, fc-correctcd from observed 0.5—2 keV. Sec text for more details. 

° Dust-corrected SFR 

Upper limits are 2.5-a" values. 



3.1. Individually detected X-ray sources 

We use the p robabilistic matches t o op tical/IR sources 
from the 1X1 II catalog (see Section I2.2[) to identify 20 
LBGs, with redshifts ranging from z — 1.5 — 5.0, that 
wer e indi vidually detected in the Chandra 4 Ms data. 
The Kill catalog lists photometric or spectroscopic red- 
shifts for all the sources (shown in column 5 in Table 
[3]). Therefore, the redshifts determined by the LBG 
dropout technique offer another measurement for com- 
parison. The spectroscopic redshifts match those of the 
dropout technique reasonably well (8 of the 12 redshifts 
agree within |Az| < 0.5). The 4 cases where the spec- 
troscopic redshift differs significantly (|Az| > 0.5) from 
the dropout redshift include 3 [/-dropouts (z ~ 3.0; 
XIDs=:563, 577, and 588) and 1 i?435-dropout (z « 
3.8; XID=546), and in all these cases the spectroscopic 
redshift quality is "secure" (> 95% confidence levels 
with multiple spectral features). With the exception of 
XID=563, all of these were detected in L05 and with 
the same dropout classification. While the discrepan- 
cies in values between the spectroscopic redshifts and 
the LBG selection seem considerably higher than the ex- 
pected contamination rate (3-10%, see Section the 
redshift failure rates in these X-ray detected LBGs are 
not representative of the LBG sample in general. In fact, 
these 4 AGN are the most UV luminous of our entire sam- 
ple, and therefore, most likely to contain high-ionization 
emission lines that are not modeled by the LBG selection. 

The highest redshift LBG sample containing detected 
sources is z = 5, with a single X-ray detected source 
(XID=371). Wh ile the photometric redshift for this 
source in the IXIII catalog is z = 3.1, a more recent de- 
termination utilizing the CANDELS photometry gives a 
more likely redshift of z = 4.651°;^^ (IXue et all 12011) . 

We find that 5 of the 20 X-ray detected LBGs ap- 
pear to be in close pairs (with other LBGs, in the same 
redshift sample, within 30 kpc; marked with asterisks 



in Table O. For comparison, 512 of the 4419 LBGs 
(^ 12%) are in close pairs. Therefore, based on our 
crude estimate. X-ray detected LBGs appear to have 
a higher probability (factor of 2 enhancement) of be- 
ing in a close pair than other LBGs. However, our 
analysis does not calculate the real space pair fraction 
{i.e., the width of the redshift bins, listed in Table [U 
correspond to large physical distances) nor takes into 
account the possibility of chance projections; therefore, 
wc stress that a more detailed statistical analysis is re- 
quired to measure accurate pair fractions and test the 
significance of close pair enhancement in X-ray detected 
LBGs. In one case, XID=386, the other LBG is sepa- 
rated by 0.8"(at z = 3, this corresponds to a physical 
separation of 6 kpc; Figure |4]). Inspecting the HST ACS 
zgso image for this source, we observe two bright peaks, 
surrounded by diffuse extended emission. Given that the 
observed sizes of LBGs at z ^ 3 are typically sinaller 
(^ 2-3 kpc; [F crguso n'etal|[20p: iTruiillo et afl [20061: 
iFranx et aL|[2008t iMosleh et al.ll2011[ ). it is unlikelv that 
the LBGs belong to a single galaxy larger than 6 kpc. 
Rather, we speculate that these two separate galaxies 
are undergoing a merger. These sources also correspond 
to a detected radio source of ~ 8a. However, none of the 
other 4 close pairs were radio detections; in these cases, 
the angular separations are 3-4", corresponding to phys- 
ical separations > 20 — 30 kpc at z = 3.0 and 3.8. While 
further study of these pairs is beyond the scope of this 
paper, we note this potentially relevant observation and 
conjecture that the enhanced star formation and AGN 
activity could be related to i nteractions in these LBGs 
(see also iKocevski et aLll2012D . 

3.2. X-rays and star-formation 

Concentrating on the X-ray properties related to star- 
formation in LBGs, we split our five z < 4 LBG sam- 
ples into two separate SFR bins for stacking: SFR/(M0 
yr^i ) = 5 - 15 (medium; see Figure [l]) and SFR/(M0 
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Fig. 6. — We compare the 2-10 keV (rest-frame) luminosity ver- 
sus SFR for different galaxy populations to show that UV-selected 
galaxies appear to dev iate towards higher Lx- Local galaxies 
(gray points, from [Colbert claL 1004) and LIRGs/ULIRGs (or - 
ange triangles fromTwas awa et al.ll2009l and lLehmer et 311120101 ). 
Stacked medium SFR LBGs (open green squares), stacked high 
SFR LBGs (open blue stars), two local analogs (cyan diamonds, 
Haroll; Grimes et al. 2007 and VV114; Grimes et al. 2006) all 
have higher X-ray luminosities per SFR compared to other galaxy 
populations. The LBG symbol sizes correspond to the redshifts, 
smallest (2 = 1.5) to largest (2 = 3.8). Solid gray line (shaded 
region) shows correlation (scatter) based on low SFR galaxies and 
solid orange line shows fi t to data, including IR-selccted galaxies, 
from ILehrner etaTIIMlOl . 

yr-1 ) = 15-100 (high). We do not include the low SFR 
galaxies (SFR =1 — 5 yr~^ ) in our analysis because 
the LBG samples suffer from incompleteness here (see 
Figure [2|) and there are no stacked detections, with up- 
per limits providing only weak constraints. Stacking of 
the high SFR LBGs all yield detections (with 2-10 keV 
X-ray luminosities ranging from 1.6-3.0 x iC'^erg s~^), 
while only the z w 1.5 and 3 samples are detected for 
the medium SFR bins. The results from this analysis are 
summarized in Table [5l Our least significant detection 
(at 2.6-ct) is the z = 1.5 sample, whose stacked 0.5-2 
keV image is shown on the left side of Figure [S] Within 
each redshift and SFR bin, the mean and median values 
of the SFRs (given in columns 6 and 7 of Table [2]) agree 
fairly well, suggesting that the X-ray stacking results are 
representative for the whole sample and not skewed by 
galaxies in a particular SFR regime. 

In Figure [6l we show the X-ray /SFR relation (Lx 
vs. SFR) for the stacked LBGs (open blue stars; sizes 
relate to redshift, smallest to largest range from z = 
1.5—3.8) as compared to othe r samples of s t ar-for ming 
galaxies: local galaxies from iColbert et al.l ()2004f ). lu- 
minous infrared and ultraluminous infr ared galaxies 
(LIRGs/UL IRGs; orange t riangl es) from iLehmer et al.l 
(|2010) and llwasawa et all (12003). VVl 14 and Haroll 
are two local LBG analogs (jGrimes et al.l 12006. 200l) 
and shown as labelled cyan diamonds in Figure [51 The 
stacked LBGs appear to have a similar X-ray/SFR ratio 
to the low SFR local galaxies (gray points). However, 
all of the UV-selected samples (the local analogs, VV114 



and Haroll, and the stacked z = 1.5 — 4 LBGs) appear 
to have slightly higher Lx per SFR. The gray line (and 
shaded region) follow the Lx versus SFR relation (and 
1-a scatter) derived in ILehmer et al.l (|2010f ) for the gray 
points, (hereafter, the local relation). The high SFR end 
of this relation is difficult to study since few local galax- 
ies have such high SFRs. However. ILehmer et al.l ()2010f ) 
find that the infrared-selected galaxies (LIRGs/ULIRGs) 
deviate from the local relation, with lower Lx per SFR 
(i.e., they are better fit by the orange solid line). They 
argue that in typical galaxies, both HMXBs and LMXBs 
provide substantial contributions to the total X-ray lu- 
minosity, while high specific SFR galaxies have X-ray 
emission dominated by HMXBs alone, thereby causing a 
lower Lx/SFR (see equation [T]) . Stud ying the X-ray/IR 
correl ation in < z < 2 galaxies, iSvmeonidis et all 
(j201lD also find that the correlation deviates from the lo- 
cal correlation at high Ljr, presumably due to relatively 
high specific SFRs. Other local (z < 0.2) samples of 
UV-selected galaxies, with high specific SFRs, can offer 
important insight at these high SFRs and further inves- 
tigation of this rela tion a t the high SFR end is ongoing 
(Basu-Zych et al. i n prepf ). 

While UV-selected galaxies (locally and at z > 1.5) 
appear similarly shifted in the X-ray/SFR relation com- 
pared to other galaxy populations, it has not been firmly 
established whether the local X-ra y/SFR relati o n can 
be applied to high redshift galaxies. iCowie et al.l ()2011l ) 
claim that the ratio of X-ray luminosity to UV luminos- 
ity does not evolve from the local relation, if you assume 
a constant dust attenuation factor (with values between 
3-5) for stacked z =1-6 LBGs. In Figure the stacked 
LBGs (open green squares and blue stars) show subtle 
variation with redshift (symbol size ranges from z = 1.5, 
smallest, to z = 4, largest). However, extending our 
analysis to z < 1.5 star-forming galaxies and taking dust 
attenuation corrections into account, we further examine 
the X-ray/SFR relation for traces of evolution. 

In Figure [3 we show the mean Lx/S FR ratio of 
our s ample and further include results from iLaird et al.l 
(Mr5, filled green diamond) for SFR=5-15 Mq yr'^, 
Xaird et al.] (|2006l . filled cyan d i amon d) for SFR=15- 
100 M0 yr-i, and iReddv et al.l ()2006L open purple di- 
amonds) for SFR> 30 Mq yr'^. For comparison with 
lower redshift samples, we include the local star-forming 
galaxies from ILehmer et all ('2010') and the < z < 1.4 
star-forming galaxies from jLOS , . We not e that the high 
SFR point at z = from ILehmer et al.l ()2010t ) is dom- 
inated by IR-selected galaxies (as discussed previously, 
these appear to have lower Lx per SFR). 

Combining the information from the data shown in 
Figure [71 including our LBGs sampled by redshift and 
SFR bins (given in Table [5]) and other published data, 
we parameterize the X-ray luminosity in terms of redshift 
and SFR. Using error-weighted least-squares fitting to 
the stacked detections, we derive the following best fit to 
the available data: 



log Lx = A log(l + z) + B log SFR + C (6) 
A = 0.93 ±0.07 
5 = 0.65 ±0.03 
C = 39.80 ±0.03 
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Fig. 7. — We show the evolution of rest-frame hard (2—10 keV) X-ray luminosity per SFR between < 2 < 4 (the last 12 Gyr of cosmic 
history). The ie/t panel shows all the data, while the right panel separates the results into different SFR ranges to show redshift trends. 
Symbols mark bins of galaxies with different SFR ranges: gray circles have low SFRs, green squares show medium SFRs, blue stars show 
high SFR galaxies, and purple diamonds mark the highest SFR galaxies (SFR/M© yr~^> 30). Following the same scheme, larger fille d 
symbols show local (2 ~ 0; Lehmer et al. 2010) galaxies, and smaller filled symbols are star-forming galaxies bet ween < z < 1.4 ||L08I '). 
Stacking results from other star-forming galaxies are shown : 2 = 1 Balmer break galaxies (filled green dia mond: .Laird et al.l2005l) , 2 =1.5-3 
BX/BM galaxies (open purple diamonds: IReddv'fc Steide]||200l) . and 2 = 3 LBGs (filled cyan diamond: ILaird et al.ll2006f) . Left; The red 
hatched region shows the regime, log Lx = 41.0 — 43.5 erg s "-*-, inhabited by individually detected X-ray sources, assumed to be AGN. The 
gray curves show X-ray binary synthesis models from'FO for total (LMXB-|-HMXBs; solid), HMXBs (dashed) and LMXBs (dotted) for 
galaxies with SFRs >1 Mp) yr~-^. Right: The gray line and shaded region (shown in all panels) represent the local X-ray/SFR relation and 
its scatter (derived by IL08I) . The best fit m odel parameterization (solid curves) show weak redshift evolution, as described by Equation 
[6l X-ray binary synthesis models from |FT 2 (dashed curves) are not fits to the data and show excellent agreement with the data and our 
best-fit. 



The trend for increasing Lx/SFR with decreasing SFR is 
consisten t with what is observ ed for local galaxies. For 
example, iLehmer et al.l (|2010f ) find Lx varies with SFR 
in a similar way (within uncertainties, referring to coefh- 
cient B in their Table 4) for their local sample of LIRGs. 

In addition, we find evidence of weak redshift evolution 
(shown as solid curve in right panel of Figure [71 where 
SFR in the equation is set to the mean SFR of the dis- 
played sample). Our result (Equation [6]) is robust even 
when we restrict our analysis to the most complete sam- 
ples (high SFR LBGs), where selection effects {i.e., flux 
limits) arc minimal. Therefore, we argue that this evolu- 
tion is driven by that of physical properties {e.g., metal- 
li city, star formation history) within the galaxies. 

IFI2I perform a large scale population synthesis study, 
using the Star Track b inary population synthesis code 
dBelczynskLeLalJlSQO^ ^2008), that models XRB popu- 
lations from the first galaxies of the Universe until to- 
day. They use as input to their modeling the Millen- 
nium II Cosmological Simulation (iBovla n-Kolchin et al.l 
20091) and the u pdated semi-analytic galaxy catalog by 
Guo et al] (j2011l ) to self-consistently account for the star 
formation history and metallicity evolution of the uni- 
verse. Their models, which are constrained by the ob- 



served X-ray properties of local galaxies, give predic- 
tions about the global scaling of emission from X-ray bi- 
nary populations with properties such as SFR and stellar 
mass, and the evolution of these relations with redshift. 

The results presented in F12 correspond to a general 
population of galaxies, including a mixture of actively 
star-forming and passive early type galaxies. In this pa- 
per, we consid er the maximum likelihood model, as re- 
ported in lF12l which we adapt to the galaxy samp le prop- 
erties of our survey. More specifically, while IfT2| used all 
the galaxies in the Millennium II simulation box in order 
to derive the star-formation history and metallicity evo- 
lution, here we only took into account the evolution of 
galaxies with SFR in four bins (1-5, 5-15, 15-100, >30 
Mq yr~^). We assume that LBGs are representative of 
the galaxy population over these redshifts {z =1-5) and 
therefore any metallicity and star formation history evo- 
lution in the LBGs will also be described by galaxi es in 
the Millenium II simulation (and therefore, in the IF12I 
models). These adapted models can be directly com- 
pared to our stacking results (see Figures [7]-[ni) • 

Noting that these X-ray binary synthesis models are 
not fits to our data, there is remarkable agreement of 
these models with the data and empirical fits. The top 
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Fig. 8. — We compa re re siduals between data and X-ray binary 
synthesis models fromU^ (top panel) and empirical fits (bottom 
panel) versus redshift. The colors and symbols are as described for 
Figure[7] Top panel: Fl^ models agree excellently with the data at 
all rcdshifts but best for the low SFR and high SFR samples at 2 = 
and the SFR > 30 Mq yr"^ at z > 2.5. Bottom panel: We display 
the residuals of the empirical fit ( log L x is fit by the equation stated 
here) to the data (symbols) and [FIJ models (dashed curves; with 
the colors corresponding to the SFR bins as in the previous Figure: 
gray, green, blue and purple represent SFR/Mq yr~^ = 1—5, 5—15, 
15 — 100 and >30, respectively). 

panel of Figure |8] shows the residuals between the IF12I 
models and data versus redshift. Residuals between the 
empirical fit and the data are shown in the bottom panel 
of Figure [8] The dashed curves show t he d ifferences be- 
tween this best-fit relationship and the lF12l models, with 
the colors corresponding to SFR bins as in the previous 
Figure: gray, green, blue and purple represent SFR/(M0 
yr~^ ) =1-5, 5-15 , 15-100 and >30, respectively. 

We use the lF12l models to interpret the evolution in the 
X-ray/SFR relation. The evolution of 2-10 keV emission 
per SFR from XRBs involves several competing effects: 
star formation history, the evolution of XRB populations, 
and metallicity evolution. The contributions of HMXBs 
and LMXBs shift with the global star formation history 
of the Universe; HMXBs trace the young stellar popula- 
tions and are expected to scale with SFR, while LMXBs 
track older stellar populations which dominate once the 
young stars are extinct. This effect can be seen in the 
shape of the LMXB curve (dotted gray) in Figure[7][Ze/i(], 
which includes the global star formation history. Metal- 
licity affects the number of HMXBs and their X-ray lu- 
minosities; high mass stars, which eventually evolve into 
black hole HMXBs (more X-ray luminous than neutron 
star HMXBs) are easier to produce in high metal licity 
environments because of weaker stellar winds ()F12[ ). At 
higher redshifts, metallicities are lower and the contribu- 
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Fig. 9. — Top: The observed evolution of the normal-galaxy 
XLD, px) from our stacked LBGs (solid black points) and z < 1.4 
late-type galax y X-ray luminosity function (open gray circles; 
IPtak et al.l 120 07) shows excellent agreement with XRB models 
(solid black, dashed blue and red dotted lines for total, HMXBs 
and LMXBs, respectively; for SFR>1Mq yr~^ galaxies). Bot- 
tom: Open circles show total px (including individually detected 
sources) from stacking, which are primarily dominated by X-ray 
emission from AGN. Results from other stacking analyses are 
shown in gray symbols, as described in the legend. The curves 
show X-ray luminosity function models, as labelled, where thinner 
lines show extrapolated values. The axis on the right-side provides 
black hole mass accretion rates (see text for details) . Since no indi- 
vidually detected LBGs were present in the z > d sample, the solid 
(top panel) and open (bottom panel) black circles are identical at 
these redshifts. 

tion of HMXBs to the X-ray luminosity is higher for a 
given SFR, apparent from the gray dashed line (HMXB) 
in Figure [7] [left]. The combination of these effects on 
the X-ray luminosity per SFR is shown as the solid gray 
line in left panel of Figure [T] With additional infor- 
mation {e.g., accurate stellar masses, metallicities and 
star formation histories), we could further explore how 
well the X-ray binary synthesis models describe these 
separate effects in observed galaxies. However, in the 
absence of such data, we interpret the evolution of X- 
ray luminosity per SFR to be driven by the metallicity 
evolution of HMXBs. Given their high sSFRs (IQ-^ '^- 
lO"*", consistent with containing young stellar popula- 
tions) , HMXBs are likely to dominate the X-ray emission 
in LBGs; metallicity evolution in HMXBs causes an in- 
crease in the Lx/SFR with redshift. 

3.3. Evolution of X-ray luminosity density as related to 
X-ray binaries 

To determine the total X-ray contribution of X-ray 
binaries in the Universe, we stack the full samples of 
LBGs (including all SFRs, but excluding X-ray detected 
sources). The top panel of Figure |9] provides our mea- 
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surements of the 2-10 keV X-ray luminosity density 
(XLD), px, due to star- formation from LBGs. px was 
calculated in the following way: 

/Ox = ?ilbg(Lx) (7) 

/OO 
0(Luv)dLuv (8) 

by multiplying the average (stacked) rest-frame 2-10 keV 
X-ray luminosity, (Lx), by the LBG number density 
(by integrating the LBG luminosity function, 0(Luv); 
down to the faintest observed UV magnitude limit for 
the sample; refer to papers listed in Table [1] for lumi- 
nosity function fitting parameters). We also include px 
values determined from the X-ray luminosity function, 
fit by a Schechter function, of lat e-type, low redshift 
(z < 1.4), galaxies (|Ptak et al.ll2007l converted from 0.5- 
2 keV band into 2-10 keV using a power law SED with 
r = 2), shown as gray open circles. We compare the 
measured values to those predicted by IF12I based on to- 
tal X-ray binary contribution (black) and separated into 
contributions from HMXBs (dashed blue) and LMXBs 
(dotted red). In general, there is good agreement of the 
total X-ray binary curve with the data; the models in- 
clude all galaxies with SFRs > 1 Mq yi~^, while the 
minimum SFR varies slightly between redshift samples 
in the data (see Figure[2]). For example, at z 1.5, the 
measured value appears slightly lower than the expected 
Px from HMXBs + LMXBs, which may be caused by in- 
completeness in the z « 1.5 LBG sample. Another expla- 
nation is that z — 1.5 LBGs may not be representative of 
all high SFR (> 1 Mq yr~^ galaxies; i.e., the LBG selec- 
tion would miss dusty star-forming galaxies like LIRGs 
or low sSFR galaxies dominated by older stellar popu- 
lations). The z « 2.5 measured value appears slightly 
higher than the model, and may include contamination 
from low luminosity AGN. 

3.4. z 3 — 8 constraints on the evolution of 
supermassive black holes 

At z > 5, we are unable to obtain detections by 
stacking only sources that are individually X-ray unde- 
tected. However, by including individually detected X- 
ray sources, we can provide constraints on black hole 
growth in these galaxies. Table [6] summarizes the re- 
sults from our z > 5 stacking. The top section of Ta- 
ble [6] excludes individually detected X-ray sources, while 
the bottom section lists the results from stacking all the 
sources. We note that including X-ray detected sources 
in our stacking analysis at z > 5 only affects the mea- 
surement at 2 sa 5 (slightly, because of one X-ray de- 
tected source). Therefore, while the results at z > 6 
remain unchanged, we shift our focus from studying ac- 
cretion in X-ray binaries to that in supermassive black 
holes (SMBHs). 

In the bottom panel of Figure IHl we show our mea- 
surements of the 2-10 keV XLD, px, from LBGs. px was 
calculated by multiplying the average (stacked) X-ray lu- 
minosity, (Lx), by the number density (as described in 
Section [321 for the top panel of this figure). Here, the X- 
ray luminosities have been calculated assuming F = 1.8, 
which is appropriate for AGN (whereas, all of our pre- 
vious calculations have been using F = 1.7, see Section 
12. 3p . This observed quantity, px, can be used to estimate 



a black hole mass accretion rate density (pbh, shown in 
log units on axis at right side) for the Universe by assum- 
ing a black hole radiation efficiency (e) and bolometric 
correction (Cboi) as: 

PBH = '^—fiCboiPx Moyr-^Mpc-^ (9) 
ec^ 

where c is the speed of light, px is given on the axis on 
left , and we adopt the va lues e = 0.1 and Cboi =25 follow- 
ing [^eiiteLiiOn (|2010t) . In this figure, we show the total 
(including the single individually detected source, listed 
in Table [3]) stacked X-ray luminosity per Mpc~'^ as open 
black circles. To compare with other studies, we use the 
0.5-2 keV count ratef|__l from other high redshift stacking 
studies to convert consistently to X-ray luminosities as 
described in Section 12.31 these results a re shown as the 
gray u pward triangl e (fFiore et aI]|2012bD, dow nward tri- 
angle (jWillott 2011), small circle (Co wie et aT1l2011[ ) and 
square symbols (Treister et al. 2010). 

We note these other studies have more sources in their 
stack since they limit sources with off-axis angles greater 
than 9' versus our adopted 7' off-axis limit; however we 
have checked that using the 9' hmit docs not change the 
results of our stacking analysis by much ( < 10%). I n 
agre ement w i thiFiore et al.l (j2012bf ) , ICowie et"all (|2011D , 
and iwiiiotti mm . we do not obtain a detection for 
the stacked z ^ 5.9 sample (see right image in Fig- 
ure [5]) , which was previously reported in iTreister et al.l 
([2010). Rather our upper limits for the z ^ 6 and 7 
data ar e very similar to those of iFiore et al.l (|2012b') . 
WillotO (|201ll) uses an optimized aperture and weighting 
in order to push their upper limits slightly lower than 
the other^3 At z « 6, o ur upper limit is low er than 
the 5(7 detection claimed bv ITreister et al.l (|2010l though 
this detection has been called into question in other pa- 
pers; see Willott 2011; Cowie et al. 201^. The solid, 
dashed, and dott ed lines show the 2 -10 keV XLD evolu - 
tion model from lAird et al I (|2010( ). lUeda et al.l (|2003| ). 
and [Silverman et al.1 ()2008[ ). extrapolated beyond the 
measured redshifts (at z = 4, 3.0, and 6.0, respectively; 
displayed in thinner lines) assuming the same form and 
derived parameters. Our observations are most consis- 
tent with the XLD evolution of AGN with Lx > lO'*^ erg 
s~^ as predicted bv lAird et al.l (|2010( ) and les s consistent 
(based on the upper limits at z « 5 — 6) with lUeda et al.1 
1120031). Yet , we a re not able to definitively rule out the 
lUeda et al.l (|2003D models since extrapolation to higher 
redshifts (z > 3) may not be accurate, and at z = 3 
the observations are consistent with their model. Ad- 
ditionally, calculating px from AGN in LBGs alone may 
underestimate the XLD from the entire AGN population, 
present also in populations of galaxies not selected by the 
LBG selection {e.g., submillimeter galaxies). We may be 
underesti mating t he con tributions from Compton-thick 
AGN: iGiUi et al.l (poTTI ) detect a Compton-thick AGN 
at z = 4.76, suggesting that Compton-thick AGN may 
be important (though difficult to detect and study) for 

where 3a upper limits are given, we converted to 2.5cr upper 
limits to the count rates for consistency with our analysis. 

22 Since IFiore et al] l l2012bl ') and iWillotH l l2011i') do not include 
CANDELS LBGs, their upper limits are not as constrained as ours. 
But when we exclude the CANDELS data, our upper limits are the 
same as theirs. 
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TABLE 6 

Properties of Stacked 1.5 < 2; < 8 LBGs Related to X-ray Luminosity Density 



z # of sources ttot Not counts S/N log Lx^ log Lyv logPx'' 

(Ms) (erg s^-"-) (erg s"-') (erg s"-"- Mpc"'^) 



Total X-ray from X-ray binaries (No SFR cut; excluding individually detected X-ray sources^) 



1.5 


30 


120.0 


56.26 ± 14.36 


4.2 


41.06 


44.12 


38.19± 0.10 


1.9 


44 


176.0 


69.16 ± 17.69 


4.3 


41.20 


44.17 


38.45± 0.10 


2.5 


222 


888.0 


117.68 ± 32.89 


3.3 


41.01 


44.15 


38.90± 0.11 


3.0 


311 


1244.0 


226.38 ± 49.08 


5.4 


41.21 


44.50 


38.77± 0.09 


3.8 


1381 


5524.0 


< 222.48 


1.6 


< 40.82 


44.12 


< 38.82 


5.0 


701 


2804.0 


< 158.51 


0.9 


< 41.18 


44.18 


< 39.00 



Total X-ray from all sources (including individually detected X-ray sources ) 



3.0 


319 


1276.0 


3513.55± 1736.9 


82.1 


42.39 


44.52 


39.96± 0.17 


3.8 


1388 


5552.0 


667.99± 445.7 


7.5 


41.31 


44.12 


39.31± 0.22 


5.0 


702 


2808.0 


< 158.62 


2.0 


< 41.20 


44.18 


< 39.02 


5.9 


220 


880.0 


< 88.56 


1.0 


< 41.62 


44.10 


< 39.23 


6.8 


61 


244.0 


< 47.07 


-1.0 


< 41.98 


44.18 


< 39.86 


8.0 


44 


176.0 


< 40.03 


-1.8 


< 42.20 


44.12 


< 39.71 



^ Total luminosity in 2—10 kcV (rest-frame) band. 

^ Total luminosity density in 2—10 keV (rest-frame) band. 

^ Since z — 5.9, 6.8 and 8.0 have no individually detected X-ray sources, the information for these redshifts is the same as given below. 

We exclude z — 1.5, 1.9 and 2.5 samples from this analysis because the volume probed by those samples is small (see areal coverage 
map in Figure and since luminous X-ray sources are rare, any analysis regarding the total X-ray luminosity probed by LBGs in this 
redshift range is too incomplete to yield meaningful results. 



understanding black hole evolution in the early universe. 
4. SUMMARY AND FUTURE WORK 

Using the deepest X-ray survey to date, the CDF-S 
4 Ms data, and the most comprehensive LBG catalogs, 
we have stacked 1.5 < z < 8 LBGs. We have split our 
analysis into two separate redshift ranges: z < 4 to study 
the relation between SFR and X-rays produced by X-ray 
binaries, and z > 5 to place constraints on the black hole 
accretion history of the universe. Our main results are 
as follows: 

1. We e xpand the list of X-ray detections (found in the 
1X1 II catalog) associated with LBGs (see Table [3]), 
based on probabilistic matching between X-ray and 
optical/infrared catalogs. Within 7' of the Chan- 
dra pointing, we match 20 LBGs with individually 
detected X-ray sources, assumed to be AGN based 
on their X-ray luminosities (> 10**^ erg s~^). 

2. While the local X-ray/SFR relation does appear to 
apply to the high redshift LBGs, it also has sig- 
nificant scatter. We find that the population of 
UV-selected galaxies (including local LBG analogs, 
also selected by using the rest-frame UV) preferen- 
tially lie at the high end of the X-ray/SFR correla- 
tion scatter (see Figure O. The high SFR regime 
of this relation has not been well-studied for dif- 
ferent samples of galaxies, but these observations 
deviate from what is observe d for high SFR, IR- 
selected LIRGs and UL IRGs (|Lehmer et al.|[20Tol: 
iSvmeonidis et al.ll20lH) . 

3. We discover mild evolution of Lx per SFR between 
< z < 4 (see Figure [7]). We stacked the LBGs in 
two different SFR ranges: 5-15 M0 yr~^ (medium) 
and 15-100 M© yr"! (high). We obtain > 2.5f7 de- 
tections in the stacked high SFR bins, and for two 



bins (z — 1.5 and 3) in the stacked medium SFR 
samples. We compare our results for the high SFR 

PJaxies with XRB population synthesis models 
T2) and find excellent agreement (see Figure [H])- 
We find that the rest-frame 2-10 keV X-ray lumi- 
nosity relates to SFR and redshift in the following 
way: log Lx = 0.931og(l-f-z)-f-0.651ogSFR-h39.8. 

4. We calculate the total (including all SFRs, but ex- 
cluding individually detected X-ray sources) XLD 
for different redshift samples to co mpare with XRB 
population synthesis models ()F12| . see top panel of 
Figurelll). We find that HMXBs and total (HMXBs 
-|-LMXBs) are both consistent with the data. 

5. Stacking z > 5 LBGs did not provide any detec- 
tions (see bottom panel of Figure [H]). However, 
we use the upper limits to derive the 2-10 keV 
XLD and we co mpare with results from other stud- 
ies and models llAird et al.. 2010: Silverma n et al.l 
feOOB; Ucda et al. 200^. Our result s are s imilar 
to .Fiore et al (,2012bl . iCowie et all ()2011[ ). and 
iWillottI (I201ll): we do not dete ct LBGs at z = 6 
fcontrarv to lTreister et alll2010D . Our observations 
are co nsis tent with the models from [Aird et aQ 
()2010f ) and lSilverman et al.l (|2008f l. extrapolated to 
z = 8; our upper l imit at z ^ 5 lies below the ex- 
trapolated (z > 3) lUeda et all (|2003D model. How- 
ever, we note that our comparisons with XLD evo- 
lution models are subject to the following caveats: 
extrapolation of the XLDs beyond the redshifts 
for which they were measured may not be accu- 
rate; studying AGN in LBGs may not account for 
the complete AGN population at these redshifts 
(e.i?., AGN in submillimeter galaxies); we may be 
underestimating the contribution from Compton- 
thick AGN. 
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Future missions like the James Webb Space Telescope 
(JWST) aim to discover the first galaxies, when the Uni- 
verse was only a few 100 million years old, searching for 
very high redshift {z > 10) LBGs. Unfortunately, study- 
ing such distant LBGs in the X-rays would be impossible 
with current X-ray telescopes, and X-ray instruments ca- 
pable of such ambitious ventures exist only the imagina- 
tion, as yet. However, the future does offer other oppor- 
t unities for studying X-ray emission from distant LBGs. 

iFiore et all (|2012a| ) describe the future prospects for 
studying AGN in LBGs using future X-ray facilities like 
Athena and Super-Chandra. As for studying XRB pop- 
ulations in these galaxies, deeper X-ray observations 
would add stacked detections (where we now have up- 
per limits) and better constrain the results presented in 
this paper. We estimate, using our Equation [6l that a 
5 Ms Chandra exposure {i.e., adding 1 Ms to the current 
CDF-S observations) would provide stacked detections 
(> 3a) for all the z ~ 2.5 — 3.8 medium SFR sample, 
and 6-8 Ms Chandra exposure would provide a stacked 
detection for the 100 highest SFR LBGs at z « 5. Deep 
IR observations from JWST and Hershel would provide 
additional rest-frame optical and IR photometry to pro- 
vide better redshift determinations for the current LBGs 
at z = 5 — 8, improving the purity of the LBG selec- 
tion and offering more accurate measurements of stellar 
masses, dust attenuations, star formation histories, and 
SFRs for individual LBGs. X-ray stacking by these prop- 
erties would provide insight about how the evolution of 
fundamental properties with redshift relate to the evolu- 
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